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Abstract 

We present a review of the electronic compressibility of monolayer and bi- 
layer graphene. We focus on describing theoretical calculations of the effects 
of electron-electron interactions and various types of disorder, and also give 
a summary of current experiments and describe which aspects of theory they 
support. We also include a full analysis of all commonly- used contributions 
to the tight-binding Hamiltonian of bilayer graphene and their effects on the 
compressibility. 
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1. Introduction 

The compressibility of the electron liquid, k, is a fundamental physical quan- 
tity which can yield information about the interactions within the fluid and the 
effect of extrinsic influences such as disorder. In the clean, non-interacting limit, 
the compressibility is straightforwardly expressed in theory in terms of quantities 
which can be easily derived from the band structure. Also, the compressibility 
is experimentally accessible since it is related to the quantum capacitance of the 
electron liquid, Cq, and can be measured directly by single electron transistor 
mounted on a scanning probe microscope. Both the compressibility and the 
quantum capacitance may be computed from dfj,/dn, where is the chemical 
potential and n is the carrier density relative to charge neutrality, since 

i=n'^ and Cq = Ae'^ (1) 
K an dfi 

where A is the sample area. Therefore the study of this quantity M = dfi/dn 
is of high importance in condensed matter physics and is the principal subject 
which we will address in this review. 

Monolayer graphene and its bilayer have become the center of an intense 
amount of research because of their novel properties, the potential for advances 
in the fundamental understanding of massless and massive chiral electrons, and 
for the many ways in which these materials might be applied in the creation 
of devices. In this article, we review the current state of knowledge of the 
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compressibility and related quantites in these materials. We focus on theoretical 
descriptions of the effect of disorder, electron-electron interactions, and the 
fundamental properties of the electron liquid, but we also discuss the current 
experimental data and describe which parts of the theory these measurements 
support. The range of features is more rich in bilayer graphene so we shall 
inevitably devote more attention to this topic, but this is not to diminish from 
the importance of understanding the situation in monolayer graphene. 

To accomplish this, in the remainder of this introduction, we describe the 
pertinent single particle properties of monolayer graphene and its AB-stacked 
bilayer, including the derivation of M in the clean, non-interacting limit. Then 
in Sec. [2] we review the role of electron-electron interactions and in Sec. [3] we 
describe the modifications due to various types of disorder. In Sec.|4]we compare 
these various theories to current experimental data and finally, in Sec. [S] we 
summarize the results and discuss their implications. 

We use the tight-binding formalism to describe the properties of electrons 
in graphene. Many thorough introductions to this theory already exist P, 
but we shall describe in detail the features which we shall utilize in this review 
to ensure that all notations and conventions are properly defined. In the case 
of the monolayer, there are two inequivalent lattice sites in the unit cell, which 
we label A and B sites as shown in Fig. [TJa) . The lattice constant is denoted a 
and the nearest- neighbor (NN) distance is a/\fi. The electron hopping process 
between these two sublattices is characterized by the energy 70 which yields 
the monolayer Fermi velocity vp = \/3a'-fo / (2h) « lO^ms"^. We exclude next- 
nearest neighbor hops because they do not produce any substantial change to 
the band structure, and we do not include any sublattice asymmetry since there 
is no well-controlled way of implementing this in an experimental context. 

The Hamiltonian of monolayer graphene in one valley may be written in 
leading order in momentum as 

where the sublattice basis is {A, B} in the K valley (with ^ = 1) and {B, A} in 
the K' valley (with ^ = —1). The operator ir = -\- iky is the linear expansion 
of the transition matrix elements near the K points. The energy spectrum 
associated with this Hamiltonian is 

E^f, =vhvFk (3) 

where v = ±1 denotes the conduction and valence bands and k = |k|. 

The crystal structure of bilayer graphene consists of two monolayer lattices 
stacked such that the A sublattice of the upper layer lies directly above the B 
sublattice of the lower layer, seperated by a distance a^. This is called Bernal, 
or AB stacking and is illustrated in Fig. [TJb). A dimer bond is formed between 
these two lattices and is characterized by the energy 71 in the tight-binding 
Hamiltonian. In the bilayer, the sublattice asymmetries and next-nearest inter- 
layer hops do have a significant effect on the band structure and we describe the 
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Figure 1: The lattice structure of (a) monolayer and (b) bilayer graphene annotated with the 
various hopping elements which are included in the tight-binding Hamiltonian. The dashed 
linos and larger unfilled circles denote the lower layer, the solid lines and smaller filled circles 
show the upper layer. Intra-layer and onsite elements are shown in black, inter-layer hops are 
drawn in green (grey). 



most relevant ones now. Firstly the layer potential asymmetry, which we denote 
by the energy u, breaks the inversion symmetry in the out-of-plane direction 
and lifts the degeneracy of the conduction and valence bands at the K points. 
Since this potential can be applied by gating, this generates a dynamically- 
tunable band gap The hops between Ai and lattices, parameterized 
by the energy 73 break the isotropy of the band structure reducing the symmetry 
from full rotational to C3, and induce an overlap between the conduction and 
valence bands. The hops between Ai and A^ sublattices and between Bi and 
Bu sublattices parameterized by the energy 74 combine with the intra-layer 
potential asymmetry A to induce an asymmetry between the conduction and 
valence bands. 

Using these notations, the Hamiltonian is 



-in 
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71 
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where the sublattice basis is {B^, Ai, Bi, Au} with ^ = 1 in the K valley and 
{Ai, Ba, Au, Bi} with ^ = — 1 in the K' valley. The energies 73 and 74 give the 
velocities V3^4 = •\/373,4a/(2?i). Figure [5] shows the band structure of bilayer 
graphene in three approximations. Solid lines denote the full solution of Eq. Q 
with 70 = 3.09eV, 71 = 0.4eV, 73 = 0.38eV, 74 = 0.14eV, and A = 0.022eV 
where these values have been taken so that vp — lO^ms"^ and 71 is in line 
with standard values taken in the literature. The higher-order parameters were 
taken from an experimental study which extracted their values by fitting optical 
reflectivity spectra to the band structure of the full Hamiltonian f^. Dashed 
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Figure 2; The bandstructure of bilayer graphene. (a) u = showing all four bands, (b) 
u = lOOmeV showing all four bands, (c) Low-energy region for u = 0. (d) Low-energy 
region for u = lOmeV. The legend in (d) applies to all four panels, and 'NN' stands for 
nearest-neighbor. The gray line in (a) and (c) is the dispersion of monolayer graphene. 

lines indicate the spectrum when the asymmetry-inducing elements 74 and A are 
set to zero; and the dotted lines show the spectrum when the trigonal warping 
elements 73 are also set to zero. This final approximation is the most frequently 
used in the literature. Figure HJa) shows the hyperbolic shape of the bands in 
the gapless case, and the asymmetries induced by the higher-order tight binding 
parameters. Specifically, the trigonal warping terms push the bands towards 
positive fcj;, and while is is not shown in the figure, this asymmetry causes the 
Fermi surface to become somewhat triangular, reducing its full isotropy to a 
three-fold rotational symmetry. The combined effect of 74 and A is to make 
the conduction band slightly steeper and the valence band somewhat shallower. 
We also include the linear dispersion of monolayer graphene for comparison, 
shown with the gray line. Figure [^b) shows the same spectrum but with a 
gap generated by m = lOOmcV. The 'sombrero' shape near the K point is now 
noticable, as is the distortion caused by the higher order tight binding terms 
which make the gap slightly narrower at the positive A:^. minimum. Figure 
djc) shows the low-energy region for the gapless case, and the overlap of the 
conduction and valence bands caused by 73 is now apparent. Figure ^d) shows 
the same but for a small gap given by it = lOmeV. Most of the complexity of the 
low energy band structure is encapsulated in this plot, including the distortion 
of the band edge and all the asymmetries described above. 

There is a common simplification of the bilayer Hamiltonian which is often 
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used when only the very low density regime is under investigation [4]. By ex- 
panding the Green's function of the Hamiltonian in Eq. (|4]) with respect to mo- 
mentum, one can write down an effective two band Hamiltonian with quadratic 
dispersion. This effectively excludes the and Bi sublattices from considera- 
tion, since electrons residing primarily on these sublattices form the split bands 
with energy ~ 71 near the K points and which are therefore not involved in low 
energy processes. 

In order to derive exact expressions for M, we must first find an analytical 
expression for the density in terms of the chemical potential. Elementary con- 
siderations show us that the total number of electrons in the conduction band 
(or holes in the valence band) is related to the number of filled states as follows: 



where n — N/A, and the state occupancy at zero temperature is given by 
the function /(k) which takes the value 1 if the state with wave vector k is 
occupied and otherwise. The possible anisoptropy and non-trivial topology of 
the band structure implies that the Fermi wave vector may be multiply-valued, 
and is a function of the angle 0. The factors gs and 171, are for spin and valley 
degeneracies, respectively, and both take the numerical value of 2. Evaluating 
the radial part gives 



This integral equation (which depends on the chemical potential via the Fermi 
wave vectors) can now be evaluated to give an expression for in terms of n. 
This can then be differentiated with respect to n to give M . 

This process can be carried out analytically in a number of cases. For mono- 
layer graphene, the Fermi surface is always circular such that kF_ {&) — and 
kp+i^) = kp. Then, Eq. ([5]) becomes 
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The band structure also gives us ^ = hvpkp = fivp^/wn so that 




(7) 



Similarly, the quadratic band structure yields a constant value for M: 




Figure 3; M = dfi/dn for clean bilayer graphene. (a) The gapless case, also with M"^ for 
monolayer graphene. (b) u = lOmeV, and (c) u = lOOmcV. The legend in (c) shows the 
approximation for the bilayer Hamiltonian used in each case and applies to panels (a-c). (d) 
Low-density region for the full Hamiltonian and three different gap values (corresponding to 
the solid lines from the other three panels). 



The four-band model for bilayer graphene with 73 = 74 = A = also gives an 
analytical result. However, in this case, the Fermi energy may be ring-shaped 
for low density, and for high density both the low-energy and split bands are 
populated. Therefore M shows two step-like features as the Fermi energy moves 
between these regions: 



V 9u ^ 



where we use the notation g„ = jf + for brevity. The first term of this ex- 
pression which contains the (5-function is present because in the gapped case, a 
discontinuity exists in the chemical potential at the band edge: For an infinites- 
imal change in density, the chemical potential jumps from the valence band 
to the conduction band, and differentiating this step gives the S function as 
shown. Including the higher-order tight binding elements means that analytical 
solutions are no longer possible, and numerical evaluation must be conducted 
instead. 



l-g„/(2^Ag„+7^/4) 



A< w2 

<X<2-ff + 

A > 272 + 



(9) 
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Figure [3] shows M in the various cases we have described. We point out 
the important features. Firstly, in the gapless case, shown in Fig. [3lja), the 
monolayer (gray line) M is substantially larger than the bilayer and diverges as 
1/n as n — )■ 0. In the NN approximation for the bilayer, M goes smoothly to 
the finite value h^Vpir /'yi at n = 0. The descrease in M at finite density reflects 
the non-par abolicity of the band structure. When the second-order hopping 
elements are included, a sharp peak is present at n = with dips at small density 
on either side. This is due to the change in topology of the Fermi surface with 
finite 73 and the region where M decreases relative to the NN case corresponds 
to the density range where Fermi surface has split into four pockets. At high 
density, the asymmetry induced by the combination of 74 and A manifests as a 
reduction in M for the valence band and enhancement for the conduction band. 
This also reveals that the presence of absence of 73 has little effect at higher 
density. When a small gap is opened, shown in Fig. ^h) for u = lOmeV, the 
divergence in the density of states at the band edge manifests as M — > when 
n — 0. However, we illustrate the 5-function spike as an infinitesimally-wide 
divergence exactly at n = 0. The effect of the trigonal warping is now masked by 
the band gap and the asymmetry between the conduction and valence bands is 
unmodified. When the gap is wider, shown for u — lOOmeV in Fig.[3Kc), there is 
an extended range of density (approximately \n\ < lO^^cm"^) where the Fermi 
energy is in the sombrero region. For the NN approximation, M is linear, but 
for the higher-order approximations of the Hamiltonian, the non-isotropic band 
structure creates a rather complex shape for M. In Fig. ^d) we show the low- 
density region for the gapless and gapped {u = lOmeV) cases, for the NN and 
full approximations for the band structure. This shows the dip-peak structure 
in the gapless full approximation more clearly, in addition to the complexity of 
M in the gapped full approximation. 

Finally, we note that this single particle analysis has been used to com- 
pare the quantum capacitance of bulk monolayer graphene to that of graphene 
nanoribbons Q- 

2. Electron electron interactions 

We now consider the effect of direct Coulomb interactions between electrons 
and correlations on the compressibility. The first observation Q is that in the 
absence of direct electron-electron interactions, the kinetic energy part of the 
correlation (an expression of the fact that electrons are fermionic) between two 
electrons behaves qualitatively differently for monolayer and bilayer graphene. 
Specifically, the two-particle correlation cancels exactly for monolayer (leaving 
only the contribution from the trivial single particle kinetic energy part) whereas 
it remains finite for the bilayer. This is due to the differences in the chirality 
of the two electron liquids. This cancellation in the correlation implies that 
the compressibility is determined by only the single particle kinetic energy, as 
seen in experiments In contrast, the bilayer does not show this cancellation, 
implying that correlations and the contribution from the interactions are likely 
to play a more significant role in that material. 
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To take into account the effects of interactions between electrons, various 
levels of approximation have been employed. At the Hartree-Fock level in 
monolayer graphene [loj . the renormalized M retains the n~^/^ dependence 
of the non-interacting system, but the value of M is increased relative to the 
non-interacting electrons. The overall enhancement in the extrinsic case is ap- 
proximately 20% when the parameters relevant for graphene mounted on an 
Si02 substrate are used. The authors also showed that the renormalization is 
increased when the dielectric environment of the graphene is weaker, suggesting 
that the enhancement of M will be stronger in suspended graphene. At the RPA 
level with finite doping, interactions between electrons favor states with large 
chirality ll[ which leads to a supression of both charge and spin susceptibilities. 
In this approximation, the overall sign of the compressibility remains positive 
but the interactions reduce it by up to a factor of 2, depending on the strength 
of the coupling constant (which is defined by the environment of the graphene) 
and the value of an arbitrarily defined momentum cutoff. For physically realistic 
parameters, a supression of k oc 1/M of 15% is predicted. 

In contrast, the introduction of electron-electron interactions at the Hartree- 
Fock level to ungapped bilayer graphene [l^ gives a profound change to the 
predicted inverse compressibility (proportional to M) because at low carrier 
density it becomes negative and divergent. Just as in the monolayer, the con- 
tributions from the intra-band and inter-band term are opposite in sign, but in 
the case of the bilayer at low density, the inter-band terms (which are negative) 
are stronger which leads to a negative divergence as n — ^ 0. For the mono- 
layer, the (negative) inter-band terms are never strong enough to overcome the 
(positive) intra-band contribution. Also, while the two-band model provides a 
reasonable approximation for M at low densities, the authors point out that 
the full four-band model is required at intermediate and higher densities. If the 
RPA approximation to the correlation part of the exchange is included fl^] then 
the negative divergence is removed because the contribution to the total energy 
per particle given by the correlation effects cancels the decreasing contribution 
from the exchange implying that M is always positive. The compressibility k 
is reduced with respect to the non-interacting case, but this enhancement of M 
reduces with increasing carrier density. 

Finally, we mention in passing that the compressibility may be a useful tool 
for investigating the nature of the ground state of charge-neutral (intrinsic) bi- 
layer graphene. There have been many theoretical works (Refs. and [3] 
are two of the earliest) which suggest that this system is unstable against tran- 
sitions to various broken-symmetry phases. The nematic phase and anomalous 
Hall insulator are two of the strongest candidates, and may lead to distinctive 
behavior of M and k at low density. In fact, some experiments [l^ have already 
been interpreted in this manner. 
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3. Disorder 



3.1. Introduction 

The effect of disorder is known to be crucial for determining many observ- 
able properties of graphene materials, and in many cases may mask the ef- 
fects of electron-electron interactions discussed above. For example, the Dirac 
point physics can be completely masked by charge inhomogeneity and in the 
more dirty samples, the presence of the externally-tunable band gap in bilayer 



graphene can be obscured in transport measurements [17|. The exact cause of 
the charge inhomogeneity (which manifests as 'puddles' at low average carrier 
density) is still somewhat controversial, so below we briefly outline the proposed 
mechanisms and proceed to assume the existence of the charge inhomogeneity 
as a phenomenological reality and discuss its impact on the compressibility. 

For monolayer graphene, the effect of a finite scattering rate between the 
electrons and charged impurities, the effect of ripples, and a contribution from 
the correlation between electrons has been examined [isj . and good agreement 
between the RPA theory and experimental data was found at higher density. 
However, at low density the theory does not modify the single-particle behavior 
(a divergence as n — 0) in contrast to the experimental data. 

3.2. The effect of puddles 

The 'puddles' of carriers were first observed in monolayer graphene via scan- 
ning SET microscopy l9jl and subsequently in STM measurements for mono- 
layer [3] and bilayer [20| graphene. The conventional wisdom is that these 
puddles are caused by external charged impurities (either located in the sub- 
strate, at the interface between the substrate and the graphene, or on top of the 
graphene) which have an electric field associated with them. This field produces 
a spatially-fluctuating potential which the electrons experience as a shift in the 
band structure relative to the Fermi energy which in turn modifies the number 
of carriers in that locality. This idea has been used theoretically to predict the 
size and distribution of these puddles 21,22, 23| , and these studies agree qual- 



itatively with the experimental measurements for both monolayer and bilayer 
graphene. Recently, another theory has been advanced which suggests that 
the scalar and vector potentials induced by corrugations in the graphene sheet 
may also cause puddles to form. The authors use a continuum elasticity theory 
for the graphene sheet and self-consistent Kohn-Sham-Dirac theory to calculate 
the induced carrier density. They find that there is no obvious correlation be- 
tween the corrugations and the puddle landscape and state that this is due to 
the non-local relation between the corrugation-induced scalar potential and the 
tensor field due to height fluctuations. Other work [5^ has also hinted at a link 
between charge inhomogeneity and deformations of the graphene sheet since the 
rippling of monolayer graphene on a substrate gives rise to a gauge field which 
may cause low-energy Landau levels to form, accompanied by an increase in M 
at low density. 

A phenomenological theory which accounts for puddles has been developed 
[iil . [231 for incorporating the effects of disorder in calculations of M. For a 



9 



measurement technique which simultaneously samples an area large enough to 
encompass several puddles (such as a bulk capacitance probe which samples 
the whole area between the gates, or an SET tip which samples an area with 
radius ~ lOOnm) and which couples capacitively to the electronic liquid in the 
graphene, it is sensible to assume that the measured response is the average over 
the sampled area. The idea is that each area with similar density will act as an 
independent capacitor and that M will be set by that density. Hence the total of 
all these parallel regions is their average and this suggests that the areal average 
can be replaced by an average over the density distribution P. In order to model 
this distribution we write the total local density n(r) = rto + n(r) where no is 
the average density which in an experiment is set by the external gates, and 
n(r) is the spatially fluctuating part induced by the disorder potential. We also 
assume that the density distribution is described by one parameter - a measure 
of the average fluctuation which we label 5n. The average value of M, which 
we denote M is 

/oo 
M{n,u)P{n,nf),6n)dn (10) 
-oo 

Note that this procedure contains an uncontrolled approximation, that the Bloch 
states used in the derivation of M persist in the disordered context and may 
be used as a basis over which to perform the average. Since the density in- 
homogeneity breaks translational invariance, this is not strictly true. Also, for 
monolayer graphene, this procedure is not completely defined because the n~^/^ 
divergence in M at low density means that the integral is infinite. However, and 
artificial cutoff could be introduced to rectify this. The (5-function divergence 
in the gapped bilayer is integrable and hence provides no technical difficulty. 

In the case of bilayer graphene, this procedure has been used to fit experi- 
mental data with a high degree of accuracy [3, HI] despite the phenomenological 
nature of the theory. To illustrate this. Fig. |4] shows M for different disorder 
strengths, different values of the gap, and for several combinations of the tight 
binding elements 73, 74, and A. It is known both experimentally 28] and the- 
oretically jl^ that P is Gaussian for bilayer graphene and so we write 



P{n, no, 6n) — — exp 



(n - no)^ 



2Sn^ 

As representative values of the density fluctuations, we take 6n — lO^-'^cm"^ to 
stand for graphene on an Si02 substrate, and 6n = 10"'^°cm~^ for both suspended 
graphene and graphene on an /iBN substrate. According to the literature, these 
are constitute rather conservative estimates of the density inhomogeneity. In 
Fig. Hlja) we show M for strongly disordered bilayer graphene in a wide density 
range. At higher density, the averaging procedure has very little effect and M 
is almost unchanged from M. But at low density the sharp dips in M (see Fig. 
[3]) become smoothed out into a small dip in M near n = 0. When the disorder 
is much weaker, as shown in Fig. lUJb), the sequence of dips and central peak 
at n = is restored (notice the different scale on the horizontal axis in this 
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Figure 4: Effect of averaging procedure M. All lines correspond to the legend in (a) (where 
'NN' stand for nearest- neighbor), and the parameters used in each panel are shown in the 
accompanying table. 



panel). When a small gap is opened, the extra complexity in M in the low 
density regime is obscured by the wide sampling region for strongly disordered 
samples, as shown in Fig.HJc). But for smaller disorder, Fig.|4l[d), the sequence 
of dips on either side of the n — peak are discernable. In the case of the wide 
gap with large disorder. Fig. Ufe), the effect of the second-order tight-binding 
elements is obscured and the most important factor determining the shape is 
the broadened S{n) term. 

In conclusion, the effect of the puddles on M at low carrier density in bilayer 
graphene is an intricate quantitative question in which two main effects - the 
size of the layer asymmetry u, and the size of the density fluctuations - both 
play an important role. Additionally, the shape of M is qualitatively different in 
the gapless case, and quantitatively different in the presence of a gap, depending 
on the approximation taken for the tight-binding Hamiltonian. 

3.3. Spatial fluctuation in band gap 

It has been observed in experiments for graphene on single-gated Si02 [20| 
that the size of the inter-layer potential asymmetry (the parameter u in the 
tight-binding formalism, which sets the size of the band gap) may fluctuate 
in space. Specifically, these authors fitted dl/dV maps of bilayer graphene in 
a strong magnetic field to the single particle Landau level spectrum and were 
able to extract both the magnitude and sign of u. They found variation on 
the order of tens of meV, and that the orientation of the asymmetry swapped 
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Figure 5: Comparison between gap formation mechanisms in comparison to experiments. 



between electron puddles and hole puddles. The spatial fluctuation of the gap 
was implemented in this phenomenological theory [27] to model the low-density 
region of suspended bilayers by assuming that the local gap could be written 
u{r) = uq + u{r) in an analogous way to the density fluctuations with du pa- 
rameterizing the fluctuations in u. Then, a second averaging procedure can be 
performed over the distribution of u yielding 

/oc 
M[no,Sn,u)P{u,uo,6u)du (11) 
-oo 

It was found that for these suspended samples, it was not possible to distinguish 
clearly between a gap formed by spontaneous electron-electron interaction ef- 
fects (which is encapsulated in a finite uq) and a disorder-induced gap modeled 
by finite Su. This is summarized in Fig. [5] where the best fits from the theories 
for fluctuating gap (M) and homogeneous gap (M) are compared to experimen- 
tal data for two samples. Therefore, if possible, measurements of M in samples 
with even less disorder is required to unambiguously determine the magnitude 
of the intrinsic gap. 



4. Experiments 

Experimental determinations of M or the compressibility of graphene have 
been made via scanning SET microscopy [^, [3| : and bulk capacitance measure- 
ments 2^ 33 ■ In the case of monolayer graphene [9l|, the measured M is close 
to the single particle result with a Fermi velocity of 1.1 x lO^ms"^. This corre- 
sponds to the theoretical result of persistence of the qualitative single particle 
behavior with a renormalization of the Fermi velocity. The effects of disorder 
at low density are noticable as a slight rounding of the divergence as n — >■ 0. 

There have been more experiments done on bilayer graphene. Two 2^, 3^ 



extract M from bulk capacitance measurements of dual-gated bilayer sheets in 
close contact to the gates. This geometry induces a large degree of disorder from 
charged impurities and the dominant effects shown in these experiments are the 
broadening of features due to the puddle averaging described in Sec. 13.21 The 
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asymmetry between the electron and hole bands caused by 74 and A are also 
clearly seen in [2^. In contrast, scanning SET microscopy of suspended bilayer 
graphene [3| accesses much cleaner samples. It is claimed that the increase in 
M measured near n = in these experiments comes from the formation of a 
broken symmetry ground state such as a nematic phase 1^ or an anomalous Hall 
insulator [l^. However, other features in the theory described in Sec. [5] (such as 
the negative divergence of in the ungapped bilayer, or quantitative estimates 
of the size of the renormalization of M due to electron-electron interactions) have 
not been observed. 



5. Conclusions and discussion 

The compressibility of an electron liquid and its related quantities M = 
dfi/dn and the quantum capacitance are important objects of study in con- 
densed matter physics since they give information about the intrinsic nature of 
the electron liquid and its interactions with external fields. In this review, we 
have described the band structure and detailed features of dfi/dn in monolayer 
and bilayer graphene, including the effects of electron-electron interactions and 
disorder. We have also reviewed the current status of experimental work. The 
single particle picture for dfi/dn is rather clear for monolayer graphene and for 
bilayer graphene at high carrier density. But for the bilayer in the low density 
limit, there are qualitative differences between the various commonly-used ap- 
proximations to the Hamiltonian. In particular, when the next-nearest neighbor 
hoppings are included in the gapless case, significant additional structure ap- 
pears at experimentally relevant density. Comparison with current experiments 
shows that disorder (specifically, the existence of inhomogeneity in the charge 
density landscape) is the dominant source of deviations from the single particle 
picture in graphene samples in contact with a substrate, while in suspended 
samples it is possible that the interaction effects are being observed. 

We gratefully acknowledge support from US-ONR, NRI-SRC-SWAN, and 
LPS-NSA. 
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